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We derive semiclassical quantization conditions for systems with spin. To this end one has 
to define the notion of integrability for the corresponding classical system which is given by a 
combination of the translational motion and classical spin precession. We determine the geometry 
of the invariant manifolds of this product dynamics which support semiclassical solutions of the 
wave equation. The semiclassical quantization conditions contain a new term, which is of the same 
order as the Maslov correction. This term is identified as a rotation angle for a classical spin vector. 
Applied to the relativistic Kepler problem the procedure sheds some light on the amazing success 
of Sommerfeld's theory of fine structure [Ann. Phys. (Leipzig) 51 (1916) 1-94]. 
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Semiclassical methods for multi-component wave 
equations have been a topic of constant interest over the 
last decade, both for their physical applications and the 
mathematical structures behind them ||l|, |[ ||, In 
a seminal article Littlejohn and Flynn y] summarized 
some of the previous efforts in this direction, stressed 
the importance of geometric or Berry phases in this con- 
text and developed a general quantization scheme. Their 
method, however, does not cover situations in which 
the so-called principal Weyl symbol of the Hamiltonian 
has (globally) degenerate eigenvalues. But this problem 
shows up for the Dirac equation, as we will explain be- 
low. It was emphasized by Emmrich and Weinstein |^ 
that in such a situation integrability of the so-called ray 
Hamiltonians (which in our case will be given by H'^ and 
H~ defined in eq. (||) below) is not a sufficient condition 
that allows for an explicit semiclassical quantization. We 
discuss this problem for the particular case of the Dirac 
equation, but our method also translates to more general 
situations. 

The semiclassical analysis of the Dirac equation was 
started by Pauli who showed that the rapidly oscillat- 
ing phase of a WKB-like ansatz has to solve a relativistic 
Hamilton-Jacobi equation. Later Rubinow and Keller 
1^ related the amplitude of the semiclassical solution to 
classical spin precession (i.e. Thomas precession [Q). So 
far, however, all these efforts did not result in general 
semiclassical quantization conditions as they were put 
forward by Keller for the Schrodinger equation In 
this article wc present the main steps in the derivation 
of such conditions, finally leading to eq. (|l^) below. To 
this end it will be necessary to extend the notion of inte- 
grability, see , from Hamiltonian systems to a certain 
skew product flow which arises naturally in the semiclas- 
sical treatment of the Dirac equation. We also illustrate 
our method for an example, namely the relativistic Ke- 
pler problem, which yields Sommerfeld's fine structure 
formula. 



We first briefly summarize the determination of semi- 
classical wave functions for the Dirac equation. Details 
can be found in |6[ U, |l^ . Consider the stationary Dirac 
equation iJo^' = E'^ with Hamiltonian 



Hv) = COL 



-V 



A{x) + (3mc^ + e(j){x) (1) 



defined on a suitable domain in L^(M'^)cg)C"'. It describes 
the motion of a particle with mass m, charge e and spin 
i in electro-magnetic potentials and A. The Dirac 
algebra is realized by the 4x4 matrices 







cr 

(T 



and (3 ■ 



l2 

-I2 



(2) 



where cr is the vector of Pauli matrices and I2 denotes 
the 2x2 unit matrix. We make a semiclassical ansatz 
of the form 



*(a?) 



,k>0 



(3) 



with a scalar phase function S and spinor-valued ampli- 
tudes afc. Inserting this ansatz into the Dirac equation 
and sorting by orders of h in leading order one finds 



[i7D(V5,a;)-S] ao = 
with the matrix-valued function 



(4) 



H^ip, x) = ca l^p A{x)j + (3mc 4- e(j){x) , (5) 

on classical phase space. The system (^ of linear equa- 
tions only has a solution with non-trivial ao if the ex- 
pression in square brackets has an eigenvalue zero, i.e. 
if S solves one of the two Hamilton-Jacobi equations 
H^{WS, x) = E with classical Hamiltonians 



H^{p,x) ^ e(j)± \lc^ {p--^A^ +rrfid^ (6) 
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for particles with positive and negative kinetic energy, re- 
spectively. From standard Hamilton- Jacobi theory, see 
e.g. we conclude that the rapidly oscillating phase of 
the wave function can be determined by integration 
along solutions {P±{t), X ±{t)) of Hamilton's equations 
of motion generated by the Hamiltonians (^). Locally 
we have P±{t) = VS'=^(X±(t)), and thus 



S^{x) = S^iy) 



(7) 



where we denote by y = X±(0) the (arbitrarily cho- 
sen) starting point of integration. If we set ^ := P±(0) 
we can also write [P ±{t) , X ±{t)) = 0^±(^, y) with the 
Hamiltonian flows . The eigenspaces corresponding 
to the eigenvalues H'^{p,x) of Hd{p,x) have dimen- 
sion two and we denote by V±{p,x) the 4x2 matrices 
of orthonormal eigenvectors, i.e. VljKf = I2 = VlV-, 
VlV- = = VlV+ and V+vl + V-Vl = I4, see @ 
for details. For concreteness we now seek a semiclassical 
wave function corresponding to the classical dynamics 
with , and in order to simplify notation drop the in- 
dex Since eq. (^ is a matrix equation it does not 
only require S to solve the Hamilton-Jacobi equation, 
but also ao to be of the form ao{x) = V{\7S,x)b{x) 
with a C^-valued b. 

An equation for b can be derived from the next- 
to-leading order equation, obtained when inserting the 
semiclassical ansatz (^) into the Dirac equation, by mul- 
tiplication with Vl from the left, cf. §, |, |l§, 



{VpH)V^b + -crB(V,5, x)b 



(8) 



-V4VpH{y^s,x)]b = o, 



3{p,x) 



e{e + mc^) 



p--A] xE B. (9) 



Here we used the abbreviation e := V {cp — eA)^ + m'^c^, 
and E{x) — ~'SJ(j){x) and B{x) — Vx A{x) denote the 
electric and magnetic fields, respectively. Viewed as an 
equation along the orbit (/>|f(^,y), the first term in (||) 
constitutes a time derivative along the classical transla- 
tional dynamics which we shall denote by a dot. The 
solution of (H) with vanishing 3 is known to be given by 

see e.g. §, and thus the ansatz b ^ . 
leaves us with the spin transport equation 



-<TB(0*,(|,y))u = O. 



(10) 



The solution of ( p^ ) can be written as u{t) = 
d{^, y, t)u{0) with an SU(2)-matrix d{^, y, t). We explic- 
itly indicate the dependence on the initial point (^, y) of 
the classical trajectory along which we integrate until 
time t. Through the covering map : SU(2) S0(3) 
we can associate with the spin transporter d a rotation 



matrix R(^,y,t), and one easily verifies that s(t) := 
Vl t) •s(O) solves the spin precession equation 

s = B(0*,(|,y)) xs (11) 

on the two-sphere S'^ (i.e. s e M'^, \s\ — 1). This is the 
equation of Thomas precession Q thus emerging from 
a semiclassical analysis of the Dirac equation. It turns 
out that all properties of the semiclassical wave function 
~ aoexp(-^S') can be determined from the solution 
(/)^(^,y) of Hamilton's equations of motion and the so- 
lution s{t) of eq. dO). Thus the skew product 



i;*(ty,5(0)) (^^(^,y),i?(|,y,t)s(0)) 



(12) 



which defines a flow on the extended classical phase space 
^2d ^ should be considered as the classical dynamical 



system corresponding to the Dirac equation, cf. [^j |lj . 

The key question in semiclassical quantization is now 
whether it is possible to find a single-valued wave func- 
tion ^/ ag exp{j-S) which solves the above equations. 
Let us briefly recall the procedure in the spinless case 

In standard semiclassics for the Schrodinger equation 
one invokes integrability of the classical flow ^Jif : Besides 
the classical Hamiltonian H —: Ai there are d — 1 fur- 
ther conserved quantities, A2, . . . , Ad (for a system with 
d degrees of freedom; we only specialize to d = 3 later) 
with mutually vanishing Poisson brackets, {Aj,Ak} — 0. 
Then the Theorem of Liouville and Arnold, see ||, chap- 
ter 10], guarantees that a (compact and connected) in- 
variant level set {{p,x)\A = const.} has the topology 
of a d-torus T'' on which the flows 4>\_^ , ■ ■ ■ , (t>Aa gener- 
ated by Ai, . . . ,Ad commute. By integration along the 
flow lines of , . . . , - analogous to the integration 
along in (M) - this allows for a definition of the phase 
function S which is unique up to the contributions of 
non-contractible loops. Demanding single- valuedness of 
the semiclassical wave function ^ ~ aoe^i^ yields the 
Einstein-Brillouin-Keller (EBK) quantization conditions 

pdx = 2Tih [rij + y) , n.j eZ, (13) 

Cj 

where {Cj\j — l,...,o?} denotes a basis of non- 
contractible loops on the torus characterized by the ac- 
tion variables Ij = §^ pdx. The number /ij S 
{1,2,3,4} is the Maslov index, see [|l^, of the cycle Cj 
which, roughly speaking, counts the number of points 



along Cj at which the pre-factor y det ^ becomes sin- 
gular. All these terms also appear in the situation with 
non-zero spin, and we now have to examine how the spin 
contribution modifies this picture. 

When we include the spin contribution c?(^,y, i) the 
situation becomes more complicated and integrability of 

will, in general, not be a sufficient condition to al- 
low for an explicit semiclassical quantization. This can 
be seen as follows: Transporting the spinor-valued am- 
plitude u along a closed path Cj on a Liouville- Arnold 
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torus the initial and final value, Mi and Uf, respectively, 
differ not only by a phase but are related by an SU(2)- 
transformation, Uf = djUi, dj G SU(2). Mathematically 
speaking, we are considering a connection in a C^-bundle 
with SU(2)-holonomy. If there was only one such loop, 
as in a system with one translational degree of freedom, 
we could choose Ui to be an eigenvector of dj , thus reduc- 
ing the SU(2)-holonomy to a simple phase factor. How- 
ever, for d > 2 degrees of freedom this is impossible 
since the holonomy factors for different loops arc, in gen- 
eral, given by non-commuting elements of the holonomy 
group SU(2). This is a general problem in semiclassics 
for multi-component wave equations with globally de- 
generate eigenvalues of the principal symbol, as was em- 
phasized in a general setting by Emmrich and Weinstein 

i- 

In our situation of semiclassics for spinning particles 
we will solve this problem by imposing additional condi- 
tions on the "field" B, which generates the classical spin 
precession (pi]). From a physical point of view it is not 
surprising that we need a stronger condition than just 
integrability of the translational dynamics 0^; since we 
identified the skew product ( p^ ) as the classical dynam- 
ics corresponding to the Dirac equation, we should also 
say under which circumstances we want to call the spin 
dynamics (or rather the combination of translational and 
spin dynamics) integrable. We do this by the following 
definition. 

Definition The skew product is called integrable, 
if (i) the underlying Hamiltonian flow (j)^^ is integrable 
in the sense of Liouville and Arnold and (ii) the flows 
, . ■ • , can also be extended to skew products l^i* 
on X (V;*! = X,\\) with fields Bj fulfilUng 

{A„Bk} + {B,,Ak} - Bj xBk = 0. (14) 

Condition (jlj) plays the same role as the condition 
{Aj, Ak} = does in the scalar case; it guarantees that 
all skew products IJi* commute . Under these condi- 
tions we are able to prove the following theorem. 

Theorem // the skew product flow yjj is integrable, the 
combined phase space R^'' x can be decomposed into 
invariant bundles Tg — > T'' over Liouville- Arnold tori 
T"^ with fiber . The bundles can be embedded in T'^ x 
such that the fibers are characterized by the latitude with 
respect to a local direction n{p,x), i.e. 

Te ^ {(p, X, s) e T'' X I <(s, n(p, x)) ^ 9] . (15) 

The proof of this theorem will be given elsewhere |TJ] . 
The geometry of the the invariant sets Tg is illustrated in 
figure 0: a Liouville- Arnold torus is sketched as a 2-torus; 
at two different points we show the attached sphere to- 
gether with the local axes n and a corresponding parallel 
of latitude. 




Figure 1: The invariant manifolds Tg of Y^^, see 
( p^ ) , are given by tori T'' to which at each point is 
attached the set of all points on the two-sphere S'^ 
sharing a fixed latitude 9 with respect to a varying 
axis n{p, x). 

If the skew product flow yjj is integrable, the theo- 
rem allows us to construct semiclassical wave functions 
which imply generalized quantization conditions involv- 
ing the spin degree of freedom. We briefly sketch the 
construction and then state the quantization conditions. 

As in the case without spin we define the semiclassi- 
cal wave function by integration along the flow lines of 
<j)\^ , . . . , <j)\^ . In addition we choose the C^-valued part 
u such that it is an eigenvector of (Tn{p, x) at each point 
of the Liouville- Arnold torus T''. (This is only possible if 
the skew product Y*^, and not just the Hamiltonian flow 
<j)\j, is integrable .) Then the semiclassical wave func- 
tion is unique up to the contribution of non-contractible 
loops on T**. Transporting a classical spin vector along 
such a loop Cj by a combination of the (commuting) skew 
products Yc\i, . . . , Yc\^, one finds that it is rotated by an 
angle aj , while integrability of 1^* ensures that it stays 
on the same parallel of latitude. Consequently, the semi- 
classical wave function is multiplied by a phase factor 
the sign depending on whether we have chosen 
u to be an eigenvector of crn with eigenvalue -1-1 or —1. 
Demanding single-valuedness of the wave function, the 
total phase change when moving along a loop Cj has to 
be an integer multiple of 27r, yielding the quantization 
conditions 

j pdx^ 2nh (^Uj + Y + '^''^) ' f^^^) 

where in addition to the terms in ([l3| ) the spin contribu- 
tion with the spin quantmn number nig = enters. 



3 



We remark that analogous quantization conditions 
can be derived for the Pauli equation There we 

can also choose to describe particles with arbitrary spin 
s S 5N0 by replacing the Pauli matrices <t with a 
higher dimensional irreducible representation of su(2). 
This changes neither the corresponding classical system 
(which is always given by a skew product on R^'^ x 5^) 
nor the construction of the semiclassical solutions; only 
in the quantization conditions ( p^ the spin quantum 
number rris then takes the values — s,— s+l,...,s. 

We conclude by illustrating these new quantization 
conditions for a famous example, namely Sommcrfeld's 
fine structure formula [p^ . To this end we have to quan- 
tize the relativistic Kepler problem with classical Hamil- 
tonian 

, 

a;) = --p-r + Vc^P^ +™^c4. (17) 

m 

The problem can be transformed to action and angle 
variables, see e.g. and the new Hamiltonian de- 

pends only on the two action variables Ir and L. Here 
Ir denotes the action variable corresponding to a radial 
cycle (e.g. from perihelion to aphelion and back), and 
L is the modulus of angular momentum L = x x p. 
In 1916 Sommerfeld quantized this system using the old 
quantum theory, since quantum mechanics was still to be 
invented, not to think about spin or the Dirac equation. 
Accordingly, he chose the quantization conditions 

Ir — hur and L ~ hi (18) 

with integers G No and / G N. More than ten years 
later it was confirmed that the energy levels resulting 
from these conditions are exactly the same as one finds 
by solving the corresponding Dirac equation |l^, |l^. 
This is insofar surprising as the Dirac equation not only 
includes relativistic effects, but also takes into account 
spin-orbit coupling, which Sommerfeld could not know 
about. Quantizing the problem with the new conditions 
(0) yields 

I^^n^nr + -±-^ and L^h[l + -±—^ 

(19) 

with integers and / and a Maslov contribution of ^ for 
both variables. For the spin rotation angle we find 
ai, = 27r for any spherically symmetric system [ p^ . In- 
triguingly, for the relativistic Kepler problem Ur is also 
given by 27r! Therefore, the conditions (|l^) and Som- 
mcrfeld's method yield the same values for Ir and L, 
thus leading to the same energy levels. A careful anal- 
ysis of the values that rir and I can assume (one finds 
rir > and I > \^\) shows that with the semiclassical 
quantization scheme developed here one also obtains the 
correct multiplicities, which Sommerfeld was unable to 
extract with his method. 

Summarizing, we can say that, by a freak of nature, 
Sommerfeld was able to obtain the correct energy levels 



of the Dirac hydrogen atom because, roughly speaking, 
the corrections due to wave mechanics (the Maslov term 
i) and those due to the spin i of the electron cancel for 
this particular problem. 
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